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We present a systematic study on the stability and the structural and electronic properties of 
mixed molybdenum-niobium disulphides. Using density functional theory we investigate bilateral 
doping with up to 25 % of M0S2 (NbS2) by Nb (Mo) atoms, focusing on the precise arrangement of 
dopants within the host lattices. We find that over the whole range of considered concentrations, 
Nb doping of M0S2 occurs through a substitutional mechanism. For Mo in NbS2 both interstitial 
and substitutional doping can co-exist, depending upon the particular synthesis conditions. The 
analysis of the structural and electronic modifications of the perfect bulk systems due to the doping 
is presented. We show that substitutional Nb atoms introduce electron holes to the M0S2, leading 
to a semiconductor- metal transition. On the other hand, the Mo doping of Nb2, does not alter the 
metallic behavior of the initial system. The results of the present study are compared with available 
experimental data on mixed MoS2-NbS2 (bulk and nanoparticles). 



I. INTRODUCTION 

Layered transition metal dichalcogenides belong to a 
well defined chemical and structural family character- 
ized by strong covalent intralayer bonding and weak Van 
der Waals interactions between adjacent layers. 1,2 The 
layer-type structure and the existence of a wide struc- 
tural gap between planes of metal dichalcogenides offer 
the possibility of an easy intercalation with a variety of 
atoms and molecules Despite their structural similar- 
ity, metal dichalcogenides present a wide variety of elec- 
tronic behaviors going from insulators (HfS2) to semi- 
conductors (M0S2), semimetals (TCS2) and real metals 
(NbS2). 4 However, it has been shown that the specific 
electronic properties of metal dichalcogenides can be eas- 
ily modified by low doping levels. 5,6 

During the past decades, a number of experimental 
and theoretical studies have been dedicated to these ma- 
terials, mainly due to the wide prospects of their po- 
tential uses. The interest in this class of materials has 
been recently renewed through the synthesis of number 
of nanosized forms, such as fullerenes and nanotubes, 
based on d-metal disulphides. Their exceptional me- 
chanical properties and a large surface to volume ratio 
makes them good candidates for new tribological and 
catalytic applications^^ Besides the pure metal disul- 
phide fullerenes and nanotubes, various "mixed" materi- 
als and composites have been recently synthesized, such 
as Nb^Wi-^Ssr^ and Ti-doped MoS 2 nanotubes^ 
The synthesis of such mixed systems might be considered 
as an extension of solid solutions well known from the 
bulk form into the field of nanomaterials, and represents 
a promising way of creating functionalized nanoparticles. 

In the class of mixed metal dicalcogenites, the M0S2- 
NbS2 system presents new promising technological po- 
tential. Molybdenum disulphide (M0S2) is widely used 



as catalyst in fuel desulphurization processes and tri- 
bological applications ^ 14 ! 15 ! 16 and niobium disulphide 
(NbS2) has attracted attention due to its optical, mag- 
netic and superconductive properties 17 i 18 The synthe- 
sis of mixed MoS2-NbS2 with no miscibility gap has been 
reported 19 and for low level Nb-doping of M0S2 particles, 
a semiconductor-metal transition has been observed. 5 
The catalytic properties of Mo x Nbi_ x S2 have been stud- 
ied in hydrogenation and hydrodesulphurization reac- 
tions. In contrast with pure M0S2, the solid solution 
Mo^Nbi-a^ (x = 0.4) has benn shown to be insensitive 
to H2S partial pressure, thus suggesting good capabili- 
ties for the conversion of high sulphur- loaded gas oils. 20 
Recently, inorganic fullerene-like Moi_ x Nb x S2 nanopar- 
ticles have been synthesized and characterized 21 and we 
expect that the catalytic activity of the mixed phase 
might be enhanced in the nanoparticles due to a high 
surface to volume ratio. 

In the context of this revived interest in mixed 
dichalcogenide based systems, the location and local 
organization of dopant atoms within the original lay- 
ered lattice play a role of particular importance. How- 
ever, current analytical techniques can not easily detect 
the exact arrangement of atoms in mixed systems, and 
so neither doping by inter-layer intercalation nor intra- 
plane substitution can be excluded j 12 i 21 i 22 Complemen- 
tary structural information can be provided by simula- 
tions using ab-initio techniques. 

In the this paper we present a density functional the- 
ory study of the mixed MoS2-NbS2 system. We consider 
bilateral doping of up to 25 % starting from pure molyb- 
denum disulphide and niobium disulphide bulk systems. 
Bilateral implies the doping of each metal disulphide by 
the metal from the counterpart disulphide (i.e. we dope 
M0S2 with Nb atoms and vice versa, NbS2 with Mo). The 
paper focus on the arrangement of dopant atoms within 



the mixed structures and on the electronic modifications 
induced in the perfect bulk systems due to doping. 



II. DETAILS OF THE CALCULATIONS 
A. Computational method 

We performed structural optimizations and electronic 
structure calculations within the framework of the den- 
sity function theory in the local density approximation 
(DFT-LDA) as implemented in the AIMPRO code^^ 
Molybdenum, niobium and sulphur pseudopotentials 
are generated using the Hartwingster-Goedecker-Hutter 
scheme. 25 The basis sets employed consisted of s, p, and 
d gaussian orbital functions centered at the atomic sites. 
For Nb and Mo atoms we have used a set of 40 gaussians 
multiplied by spherical harmonics up to a maximum an- 
gular momentum / max = 2, for S atoms a set of 20 gaus- 
sians with a maximum angular momentum / max = 1. 

Calculations for bulk M0S2 and NbS2 have been con- 
ducted using a 10x10x10 Monkhorst-Pack k-points set to 
sample the Brillouin zone. Defective structures have been 
described using a 4x4x1 supercell. 2x2x4 Monkhorst- 
Pack k-points set have been proved to be sufficient for 
accurate Brillouin zone integration in this large supercell. 
The atomic positions for all structures were optimized si- 
multaneously with the cell parameters. 



B. Structural models 

M0S2 crystallizes in several polymorphs, among which 
the most stable is the hexagonal 2H-M0S2, belonging to 
the P63/mmc symmetry point group. The lattice struc- 
ture consists of hexagonal planes of Mo atoms lying be- 
tween two hexagonal planes of S atoms, forming a S-Mo- 
S layer. Each Mo atom is covalently bonded with its six 
first neighbor sulphur atoms, each S is bonded to three 
molybdenum atoms. The unit cell contains two alternat- 
ing S-Mo-S layers with an ABA BAB stacking along the 
c axis, see Fig. [TJ In bulk 2H-NbS2 the local structure 
of the layers is similar to that of 2H-M0S2, with the Nb 
atoms 6 fold coordinated. The only difference between 
the structures is that the NbS2 layers are shifted and 
Nb atoms lie above each other, thus the stacking order 
within the bulk structure is ACA BCB, see Fig[TJ 

The experimental lattice parameters for 2H-M0S2 and 
2H-NbS 2 metal disulphides are a = 3.16 A, c = 12.29 A, 
and a = 3.31 A, c = 11.89 A, respectively^. Molybde- 
num and niobium oxidation states in the corresponding 
sulfate bulk systems are identical and the covalent radii of 
Mo (130 pm) and Nb (132 pm) are similar, implying that 
bilateral substitution might be feasible. Furthermore, the 
interlayer spacing for both metal dichalcogenides is wide 
enough to place the interstitial atom in between the lay- 
ers. Hence, in the present work we will analyze both, 
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Figure 1: (Color online) Top and side views of the 2H-M0S2 
and 2H-NbS2 crystal structures. Different models of doping 
presented on the bottom pictures: interstitial atoms in tetra- 
hedral (a) and octahedral (b) position; substitutional atom 
(c). 



substitutional and interstitial doping modes for accom- 
modating the metal atoms within the layered disulphide 
crystals. 

For the interstitial doping we have to distinguish two 
cases: the dopant atom occupies the tetrahedral site situ- 
ated between one sulphur triangle and one sulphur atom 
in the adjacent plane (position a, Fig. [1]) or the octahe- 
dral site between two S-triangles (position b, Fig. [1]). In 
the first case only four metal-sulphur bonds are formed, 
whereas in the second case six. In the substitutional con- 
figuration one metal atom is removed from the metal 
disulphide (MS2) layer and the vacant site is occupied 
by a dopant atom, which thus forms six metal-sulphur 
bonds (position c, Fig. [1]). 



III. RESULTS AND DISCUSSION 
A. Bulk systems 

The optimized bulk cell parameters are in excellent 
agreement with the experimental values: a = 3.15 A, c 
= 12.29 A for MoS 2 and a = 3.32 A, c = 11.92 A for 
NbS 2 . The Mo-S and Nb-S bond lengths in the pure 
bulk systems are found to be 2.42 A and 2.50 A, which 
is close to the reported experimental values respectively 
of 2.41 A and 2.47 A. 4 

The electronic structure of bulk metal dichalcogenides 
have been extensively studied both experimentally and 
theoretically. 4 » i 17 i 26 In Figj2] we present the electronic 
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Figure 2: DFT-LDA electronic band structure for pure M0S2 
and NbS2- The origin of the energies is the Fermi level. 



band structures for M0S2 and NbS2 bulk species. The 
deepest band in the plot corresponds to the S 3s states 
and is separated by a wide gap from a broad valence 
band. The bottom of the valence band corresponds to 
the hybrid 3p-S and 4d-Mo(Nb) states, and the upper 
part to the 4d z 2 Mo(Nb) states. As one can notice, the 
overall structure and dispersion of the bands for the M0S2 
and NbS2 bulks are rather similar. The main difference 
between the two systems comes from the different num- 
ber of valence electrons between the Mo and Nb atoms. 
The number of valence electrons in M0S2 is enough to fill 
completely the valence band and thus M0S2 is a semicon- 
ductor with an LDA indirect gap of 1.17 eV (experimen- 
tally measured^ as 1.23 eV) and a direct gap at the T 
point of 2.02 eV (experimentally 27 1.74 eV). However, 
NbS2 has one electron less per metal atom so that the 
top of its valence band is half-filled. Consequently, NbS2 
has a metallic behavior. Our results for bulk systems 
agree well with the previously reported calculations and 
available experimental data i 28 i 29 i 30 i 31 i 32 



B. Low concentration doping: single dopant atoms 

In order to access the most energetically favorable way 
for bilateral doping within the M0S2 -NbS2 system, we 
have calculated the formation energies of the "mixed" 
structures. For both types of doping, the formation en- 
ergy can be estimated from the general equation 

^form = ^mixed ~ [M • /lMoS 2 + N • /iNbS 2 + X • /!&} (1) 

where 



J M, N 7^ and X = for substitutional doping; 
1 M or N = and I/O for interstitial doping. 

In EqQ]E m i xe d is the total energy of the mixed system; 
M and N are the numbers of M0S2 and NbS2 units in 



the mixed systems; /imoS 2 and /iNbS 2 are respectively the 
chemical potentials for M0S2 and NbS2 using as refer- 
ence the corresponding bulk systems; X is the number 
of dopant atoms inside the cell and /id is the chemical 
potential of the dopant atom. 

As a general rule, in the synthesis process the chemical 
potential for an element as a dopant should be lower than 
that in its bulk form, otherwise this element would form 
the energetically more stable bulk phase rather than the 
mixed system. Consequently, in our study we imposed 
as boundaries the chemical potentials of dopants derived 
from the corresponding pure bulk metals //d = Mbulk or 
metal disulphides /id — Mms 2 — 2-/is, where /is is sulphur 
chemical potential as derived from its bulk phase. Thus, 
we obtain the range of formation energies calculated us- 
ing both values of the dopant's chemical potentials for 
each dopant configuration. 

Formation energies for low level doping (one dopant 
per unit cell, which corresponds to ~ 3% of doping) 
and metal-sulphur bond lengths are presented in Table H 
The inclusion of Nb dopant in interstitial sites within the 
M0S2 structure is highly energetically unfavorable since 
the formation energy for the Nb in octahedral sites is 
of the order of several eV. The Nb position in tetrahe- 
dral sites (position a, see FigH]) is not stable and after 
relaxation the Nb atom moves to a neighboring octahe- 
dral site (position b). In contrast, negative formation 
energies correspond to the substitutional model (posi- 
tion c, see FigJT]), indicating that the substitution process 
is exothermic. In this configuration the Nb-S distances 
decrease by 0.05 A compared to the NbS2 bulk, adapt- 
ing to the initial parameters of the "host" M0S2 lattice. 
The Mo-S bond lengths at neighboring sites remain unaf- 
fected. In the mirror situation, when NbS2 is doped with 
Mo, both the intercalation and substitution processes are 
energetically favorable. 33 Moreover, the range of forma- 
tion energies reported in Table [J suggests that low level 
substitutional and interstitial doping will be in competi- 
tion. The Mo-S bond distances for the substitutional Mo 
dopant increase by about 0.05 A compared to the M0S2 
bulk. However, Mo insertion between the NbS2 planes is 
characterized by a decrease of 0.02 A in the Mo-S bond 



Table I: Formation energies (eV) and average metal-sulphur 
bond distances (A) in single atom doped M0S2 and NbS2- 
Formation energies are calculated using the dopant chemical 
potentials derived from metal disulphide (E/ ormi ) or pure 
metal bulk (E/ orm2 ) phases. 
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E formi 


Ey rm 2 


Mo-S 


Nb-S 


M0S2 








2.42 




NbS 2 










2.50 


NbMo 32 S 6 4 


(inter) 


8.08 


3.80 


2.42 


2.43 


NbMo 3 iS 6 4 


(subst) 


-0.21 


-0.86 


2.42 


2.45 


MoNb 3 2S 6 4 


(inter) 


-3.24 


-6.87 


2.40 


2.49 


MoNb 3 iS 6 4 


(subst) 


-5.11 


-4.46 


2.46 


2.49 
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length. 

In Fig. [3] we present the band structure and density 
of states (DOS) for low level doping of M0S2 and NbS2 
crystals. The insertion of Nb dopant atoms within the 
M0S2 preserves the main features of the DOS of the bulk 
system. Furthermore, the interstitial Nb atom adds two 
additional nearly degenerate fully occupied energy levels 
within the band gap and a half-occupied donor level at 
the bottom of the conduction band (see Fig. [3]). The sub- 
stitution of one Mo atom by one Nb atom introduces an 
electron hole within the system and thus moves the Fermi 
level down. The Nb atom acts as an acceptor impurity in 
the molybdenum disulphide. To get further insight into 
the electronic properties of Nb-substituted M0S2 , we ex- 
amine the electron density distribution map for the the 
half-occupied electronic state associated solely with the 
doping, see Fig. [H As one may observe, the impurity 
level is not centered exclusively on the Nb atom, but de- 
localized within the radius of the third molybdenum ring 
around the niobium, thus suggesting the metallic nature 
of the mixed system. 

The substitutional doping of NbS2 by Mo acts in a 
symmetric way to the previously treated case, introduc- 
ing an additional electron instead of a hole into the sys- 
tem. Thus, the Fermi level is moved upward and the 
overall electronic structure of the system is conserved. 
Similar behavior is found for the interstitial doping with 
a slightly higher Fermi level shift. 

C. Higher dopant concentration: MoS2-NbS2 alloy 

In this section we treat the case of higher doping levels 
for the system MoS2-NbS2- Two extreme cases might be 
considered for the arrangement of dopant atoms within 
the original host matrix: dopants randomly distributed 
in the crystal forming a homogeneous alloy, or cluster- 
ized, forming local domains in the host crystal. 

The tendency of dopants to clusterize can be inves- 
tigated by considering the binding energy between two 
dopant atoms as a function of their separation. We dis- 
cuss the situations when two dopants occupy close neigh- 
boring sites or, contrary to that, are situated at the far- 
thest sites within the chosen supercell. We define the 
binding energy Ebind between dopants as: 

^bind = 2 • E\y> — ^2D (2) 

where Eid and E2D are respectively the formation en- 
ergies of the systems with one and two dopants. 

In Table [III we report formation and binding energies 
for the mentioned arrangements of dopants within the 
M0S2 and NbS2 lattices. For interstitial Nb atoms in the 
M0S2 crystal we obtain a high binding energy, which at 
two neighboring octahedral sites rises to 1.17 eV. This 
value suggests a strong driving force which promotes the 
clustering of Nb interstitials. In contrast, Nb substitu- 
tional atoms weakly interact, with binding energies of 



the order of a hundredth of an eV. Considering that the 
configurational entropy is higher for a random distribu- 
tion of dopant atoms, we thus conclude that the mixed 
system will tend to form a homogeneous alloy. A similar 
tendency for random ordering is obtained for both inter- 
stitial and substitutional Mo in NbS2 , where the binding 
energies between two doping atoms drop to about ~ — 5 
eV. 

The experimental work of Hotje et al^ has shown that 
the MoS2-NbS2 alloy can be synthesized without any mis- 
cibility gap across the concentration range going from 
pure M0S2 to pure NbS2- We have thus extended our 
investigation to higher doping levels up to a maximum 
of 25%. The percentage of doping is calculated consid- 
ering all the metal atoms within the two adjacent metal 
disulphide slabs, constituting the supercell. 

In section III Bl we have seen that the in-plane lattice 
parameters of M0S2 and NbS2 differ by about 4.8 %. 
Then, in section IIII Bl we saw that the incorporation of 
dopants causes structural perturbations related to differ- 
ent metal-sulphur bond lengths. For high doping con- 
centrations these frustrations will induce a change in the 
crystallographic parameters of the host crystal which de- 
pends on the stoichiometrics of the mixed systems. The 
data, obtained by optimizing the atomic positions simul- 
taneously with the cell parameters, are presented in Fig. 
[5] and compared to the X-ray diffraction data available 
from the work of Hotje et ali& The graph shows a good 
agreement between our theoretical results and the exper- 
imental data with a deviation of around 1% across the 
range of doping levels considered. 34 

Formation energies for the mixed MoS2-NbS2 systems 
as a function of interstitial and substitutional dopant con- 
centration are presented in Fig. [6j The range of values 
represented by the shaded areas and obtained from Eq. 
[TJ is defined by the two dopant chemical potentials of the 
pure metal disulphide and metal bulk phases. The ob- 
served trends in Ef orm for low level Nb doping of M0S2 are 
preserved as well at higher dopant concentrations show- 
ing a preferential substitutional doping. These results 
are confirmed by recent EXAFS experiments. 20 

Now we turn to the symmetric situation of Mo dop- 
ing of NbS2- As we observed before, low level doping 
by Mo substitution or intercalation within the NbS2 are 
exothermic (and competing) processes. The enhance- 
ment of dopant content leads to the overlapping of the 
regions corresponding to the formation energies for the 
interstitial and substitutional doping, see Fig. [6j Hence, 
depending upon the particular synthesis conditions, one 
may have a combination of substituted and intercalated 
molybdenum atoms within the NbS2 matrix. 

The effect of low level doping of M0S2 and NbS2 crys- 
tals on their electronic properties have been discussed in 
detail in Sec. IIII Bl Here, we mention that near the Fermi 
level the electronic structure is locally sensitive not only 
to the type but also to the degree of doping. However, 
the main tendencies described above will be preserved 
and enhanced by higher dopant concentrations. Thus, 
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Figure 3: DFT-LDA electronic band structure and total density of states for substitutional and interstitial doping of M0S2 (top) 
and NbS2 (bottom) crystals. Filled areas corresponds to the pure bulk systems. Solid and open triangles indicate respectively 
the filled and half- filled electronic bands associated with the doping of M0S2. The origin of the energies in each graph is the 
Fermi level. 



Table II: Formation energies Ef or m(eV) and binding energies Ebind (eV) for interstitial and substitutional doping of M0S2 and 
NbS2- Formation energies are calculated using the dopant chemical potentials derived from the metal disulphide. 



System (interstitial) 


Eform 


Ebind 


System (substitutional) 


Eform 


Ebind 


NbMo 32 S 6 4 


8.08 




NbMo 3 iS 6 4 


-0.21 




Nb 2 Mo 3 2S 6 4 a 


14.98 


1.17 


Nb 2 Mo3oS 6 4 a 


-0.44 


0.01 


Nb 2 Mo 3 2S 6 4 b 


16.07 


0.08 


Nb 2 Mo 3 oS 6 4 b 


-0.46 


0.03 


MoNb 32 S 6 4 


-3.24 




MoNb 3 iS 6 4 


-5.11 




Mo 2 Nb32S 6 4 a 


-0.91 


-5.57 


Mo 2 Nb 3 oS64 a 


-4.97 


-5.25 


Mo 2 Nb 3 2S 6 4 b 


-1.08 


-5.39 


Mo 2 Nb3oS 6 4 b 


-4.91 


-5.31 



a Two dopant atoms are separated by the shortest possible dis- 
tance within the supercell 

6 Two dopant atoms are as far as possible from each other within 
the chosen supercell 
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further Nb substitution in M0S2 will decrease the total 
electron concentration of the system and cause a pro- 
gressive downward shift of the Fermi level. However, the 
substitutional doping with Mo of NbS2, will increase the 
total electron concentration and lead to gradual filling of 
the d z 2 Nb band. 

Hence, over the whole range of considered concentra- 
tion, the electronic behavior of the stable mixed systems 
will be metallic. We note that our conclusions agree well 
with the experimental data. According to electrical re- 
sistivity measurements, M0S2 with 5 % of substitutional 
Nb atoms has an electrical resistivity level similar to that 
of graphite. 5 Likewise, the mixed Nb x Moi_ a: S2 nanopar- 
ticles are reported to be metallic^ 
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IV. CONCLUSIONS 

We have presented a systematic ab initio study on the 
stability, structural and electronic properties of mixed 
molybdenum niobium disulphides. Bilateral doping ef- 
fects have been investigated up to a dopant concentration 
of 25%. 

Focusing on the specific arrangement of dopant atoms, 
we observe that over the whole range of considered con- 
centrations, substitutional doping with Nb of M0S2 will 
predominate. In addition, our calculations for Mo dop- 
ing of NbS2 show that depending on the specific syn- 
thesis conditions, both interstitial and substitutional Mo 
arrangements can co-exist. The incorporation of dopant 
atoms causes structural perturbations and changes in the 
crystallographic parameters of the host crystal, which de- 
pendents on the stoichiometrics of the mixed systems. 

The difference in the number of valence electrons be- 
tween Mo and Nb atoms means that Nb substitutional 
doping of M0S2 introduces electron holes into the sys- 
tem. Mo substitutional doping of NbS2 adds electrons to 
the system, leading to an upward shift of the Fermi level. 




Figure 4: (Color online) Electron density distribution of the 
half-occupied electronic state associated with the introduction 
of one Nb substitutional atom in the M0S2 lattice. The Nb 
atom is blue, Mo atoms are red and S atoms are yellow. 



Figure 5: In-plane cell parameter for the mixed Moi-ccNb x S2 
systems as a function of Nb concentration inside the cell. Solid 
diamonds correspond to reported experimental data.— 




10 20 30 ' 70 80 90 100 
Nb concentration (%) 

Figure 6: (Color online) Formation energies for the mixed 
Moi_ cc Nb cc S2 systems as a function of the dopant chemical 
potential and dopant concentration within the cell. Open 
symbols denote the dopant chemical potential derived from 
pure metal disulphide, solid symbols stand for potential ob- 
tained from metal bulk phases. 



According to our results, these mixed disulphides have a 
metallic behavior throughout the range of stoichiometry 
considered. This result is in agreement with experimen- 
tal data available for low- level Nb doping of M0S2. 

We should note that the semiconductor to metal transi- 
tion due to Nb doping might not induce any degradation 
in the mechanical properties of the M0S2. Indeed, the 
dopant atoms are located within the planes and no inter 
plane bonding bridges are formed. Thus, mixed NbS2- 
M0S2 might present interesting capabilities for new tri- 
bological applications. 

The conclusions of the present work can be generalized 
to other mixed d metal chalcogenides presenting similar 
structural and electronic properties. 
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